December 1972 


B72- 10667 



NASA TECH BRIEF 

Marshall Space Flight Center 



NASA Tech Briefs announce new technology derived from the U.S. space program. They are issued to encourage 
commercial application. Tech Briefs are available on a subscription basis from the National Technical Information 
Service, Springfield, Virginia 22151. Requests for individual copies or questions relating to the Tech Brief program may 
be directed to the Technology Utilization Office, NASA, Code KT, Washington, D.C. 20546. 


Numerical Solution of Potential Flow Problems 
in Terms of Flux Components 


Flux components in vapor flow are numerically deter- 
mined by solution of a potential flow problem. In a two- 
dimensional Cartesian coordinate system, this flow is 
described by Laplace’s equation: 

&T 3 2 T 

+ = 0 (1) 

where T is the potential and y and z are coordinates. 
The flux components are then defined in terms of the 
potential by 
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where v and w are the flux components in the y and z 
direction, respectively, and k is the transport coefficient. 
There are several excellent numerical methods of solution 
of equation (1); however, this solution yields values of 
the velocity potential, and the velocity components must 
be obtained from these values by numerical differen- 
tiation. Consequently, it is desirable to develop a 
numerical solution for this problem in terms of the 
velocity components. 



Figure 1. 


Flux components can be determined by solving the 
differential equations which define them. These are the 
continuity principle 
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and an irrotationality condition 
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Equation (4) is equivalent to equation (1), and equation 
(5) is obtained by equating the cross-partial derivatives 
of (2) and (3). Both equations contain only first deri- 
vatives, and therefore, a centered difference approach 
is the appropriate one. 

The method of solution must be somewhat different, 
however, as the potential flow problem has split bound- 
ary conditions in both directions. Typical boundary 
conditions define either v or w at all boundaries. The 
region in which these equations apply is divided in both 
directions by grid lines. The first approach to this prob- 
lem is to determine both v and w at all intersections in 
the grid. A typical section of the grid is shown in Figure 1 . 
Finite difference analogs for the continuity equation 
and the condition for irrotational flow equation are 
centered about the center of the square unit (point 
marked by the cross), and each finite difference equa- 
tion contains four values of v and four of w. The analogs 
obtained are second-order correct. Several iterative 
methods have been developed for the solution of the 
resulting equations. This method provides an alternate 
numerical solution of potential flow problems in two 
dimensions. However, it has no particular advantages over 
the methods for solving in terms of the potential. 

A simple innovation in the method of defining the 
grid leads to a number of advantages which makes the 
method very promising. The location for the points of 
this checkerboard method are shown in Figure 2. The 
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dependent variable v is determined as the points desig- 
nated by the squares (□) (Figure 2), and w is determined 
as those designated by the triangles (A). The analog to 
the continuity equation is centered about the point 
marked by the open circle (o); and that to the irrota- 
tionality condition equation, by the filled circle (•). 
As a result, each finite difference equation contains only 
two values of v and two of w. These equations are much 
simpler, and the equations from adjacent rows can be 
combined so that all variables can be eliminated except 
the boundary conditions and values of v and w from 
the center rows of points. The matrix of coefficients 
for these equations is of the band form and can be 



solved by an algorithm. The solution has negligible 
round-off error for as many as twenty rows of each 
variable. In this manner, a direct solution is obtained 
for the potential flow problem with no iteration. 
Furthermore, even though the checkerboard method is 
also second-order correct, it has been found to have 
much less truncation error than the conventional method 
on several problems for which the exact solution is 
known. 

Notes: 

1. Information concerning this innovation may be useful 
in many areas of physics, e.g., conduction, heat 
transfer, etc. 

2. Requests for further information may be directed to: 

Technology Utilization Officer 
Marshall Space Flight Center 
Code A&PS-TU 

Marshall Space Flight Center, Alabama 35812 
Reference: B72- 10667 
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Figure 2. 


Source: D. U. von Rosenberg of 
Tulane University 
under contract to 
Marshall Space Flight Center 
(M FS-21751) 
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